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Abstract 

We give a principled approach for the selection of a boundary integral, 
retarded potential representation for the solution of scattering problems for 
the wave equation in an exterior domain. 


1 Introduction 

Let D be a bounded domain with smooth boundary, T, and set Vl = D'^ . The 
determination of the outgoing part of the solution to the classical acoustic scattering 
problem for a soft scatterer requires solving the mixed Cauchy problem for the 
wave equation: 

d^u = (?/\u in U X [0, oo), 

with tt(0, x) = dtu{Q, x) = 0 and (1) 

u{x,t) = g{x,t) for {x,t) € dVl x [0, oo). 
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Various other types of scatterers are modeled by using boundary conditions of the 
form I3{x)u{x, t) + dnu{x, t) = g{x,t), for {x, t) G dQ x [0, oo). 

Many numerical approaches to solving this problem express the solution in 
terms of the retarded potentials derived from the fundamental solution to the wave 
equation. By analogy with the time harmonic case, we define the single and double 
layer potentials by the formulae 


Sfj,{x, t) 


f \x-y\) 

J 4n\x-y\ 
r 



ny-{x- y) \ y{y,t- |a; - ^|) 
47r|x — y\ [ \x — yp 


/i(y,t - \x - y\) 
|a; - y\ 


dSy^ 


( 2 ) 


for X € Here, /i(y, s) = dsy{y, s). 

Suppose now that we represent the solution in the form 


u{x, t) = T>y{x, t) + S{afi + bfi){x, t) for (x, t) € x [0, oo). (3) 


To enforce the boundary condition on T x [0, oo), we take the limit of this repre¬ 
sentation from the exterior of D, yielding the integral equation 


yix,t) 


+ Nn{x, t) + K{a{y)y + b{y)y){x, t) 


= 9ix,t), 


(4) 


for (x, f) G r X [0, oo). Here N is the weakly singular principal part of the double 
layer restricted to T X r, and K is the single layer on T x T. 

There is a substantial literature on time-domain integral equations for scatter¬ 
ing problems, which we do not seek to review here. (See, for example, 131IH [TTl 
[5] |6l [TQ and the references therein.) We simply note that, while prior work has 
considered the use of the single layer alone, the double layer alone, as well as 
linear combinations of the two, there is relatively little discussion of a principled 
analytic approach for selecting one representation over another. Using classical re¬ 
sults from scattering theory of Lax, Morawetz, and Phillips [8], we provide criteria 
for choosing the coefficients (a(y) and b{y)) above, at least when T is a non¬ 
trapping obstacle. This is a hypothesis about the behavior of reflected light rays in 
the exterior of D, which is always satisfied if D is convex or sfar-shaped. These 
resulfs imply fhaf if fhe disfurbance g{y,t) has compacf supporf in time, fhen fhe 
solufion to the mixed Cauchy problem in equation ([T]l will eventually exhibit local 
exponential decay. 

The principle that guides our selection of representation is that, ideally, the 
source g{y, t) should display the same rate of exponential decay in time as the 
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solution itself. If this is not the case, then the mapping from fj, i-> u{x, t) in (l3]l 
must, after a certain time, involve some catastrophic cancellation, and hence can be 
expected to lose relative accuracy as time progresses. This goal can be achieved ex¬ 
actly for round spheres. For more general non-trapping obstacles, we show (mod¬ 
ulo a technical hypothesis) that certain choices of (a, b) lead to sources that decay 
exponentially, though perhaps not at the same rate as u itself. As is shown in Sec¬ 
tion m in the smooth, strictly convex case, microlocal analysis indicates that the 
choice a = 1 should be be optimal. 

Our analysis relies upon taking the Fourier transform in the time variable and 
using recent estimates for the kernel of the CFIE done by Chandler-Wilde and 
Monk, and Baskin, Spence, and Wunsch, see 11113. The results of |T| and our 
own also rely on several modern improvements to the classical results in scattering 
theory alluded to above. 

2 Time-harmonic Analysis 

Recall that 



(5) 


is the outgoing fundamental solution for (A -|- fc^), with the corresponding single 
and double layer potentials defined by: 



( 6 ) 


an 


OD 


With respect to the complex linear pairing, {u, v) = uvdS, the single layer Sk 
is self-dual and 


(V) 

dD 

If we take the Fourier transform in time: 



OO 



0 

then equation (01) becomes 
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As before the limit is taken from the exterior of D. Here and in the sequel we 
assume that g{x, f) = 0 for f < 0 ; by eausality g{x, t) ean also be taken to vanish 
for t < 0 . 

The eombined field operators, aeting on funetions on T, given by 

A(a, b;k) = ^ + Dk- iSk[ka{y) + ib{y)], (10) 

are well known to be Fredholm operators of seeond kind, provided that F is at 
least C^. For funetions a(y) > 0 and b(y) > 0 for all y € F, it ean be shown, by a 
simple integration by parts argument, that A{a, b; k) is invertible for k in the elosed 
upper half plane. To prove this we use that the dual of A{a, b; k) with respeet to 
the eomplex linear pairing is 

A'{a, b;k) = ^ + D'f,- i{a{y)k + ib{y))Sk. (11) 

Suppose that for some Im k > 0 , there is a non-trivial solution to A(a, b;k)f = 
0. By the Fredholm theory this means that there is also a non-trivial solution h to 
A'{a, b; k)h = 0. Let 

u{x) = Skh{x). (12) 

We will have oeeasion to eonsider the funetion u{x) in both the interior of D and 
its exterior. To make this distinetion elear, we will refer to the eorresponding re- 
strietion of u by u- and tt+, respeetively. 

Observe now that n_ 7 ^ 0. If it were, then |'a£)= 0 as well. In this ease 
would be an solution to (A + /c^)u+ = 0 with vanishing Diriehlet data, satisfying 
the Sommerfeld radiation eondition and it would therefore also be zero, whieh is 
impossible. Thus, when Im fe> 0 , this gives a non-trivial solution in D to 

(A -h k^)u- = 0 with dnU-(x) — i(ka(x) -j- ib(x))u- (x) = 0. (13) 

The boundary eondition is equivalent to the faet that A'{a, b; k)h = 0. 

Using the boundary eondition, and integration by parts we obtain: 


0 = J {Au-+ k"^u-)u-dx = J [k"^\u-\ —\Vu-f]dx + J dnU-U- 

dD 


D 


D 


= j[k^\u-^ — \Vu-\^]dx + i j{ka{x) + ib{x))\u-\^dS. 
D dD 


(14) 
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Let k = ki + ik 2 , and take the real and imaginary parts to see that this implies 



( 15 ) 


dD 


D 


The first relation shows that if fci 0 and /c 2 > 0, then n_ = 0 as well. If = 0 
and k2 > 0, then the second equation shows that U- = 0. If k2 = 0, but ki ^ 0, 
then the first relation implies that u_ \bD= 0; the boundary condition in ([T3] ) then 
implies that dnU- = 0 as well. It then follows from Green’s formula that tt_ = 0. 
This proves the following basic theorem: 

Theorem 1. If a{x) and b{x) are positive functions defined on dD, a compact 
-surface, then the operators A{a, 6; k) and A'[a, 6; k) are invertible for k with 


Im A: > 0. 


Remark 1. If 6 = 0, then A{a, 5; 0) has a non-trivial nullspace. Such a nullspace 
would generically destroy any possibility for exponential decay in the source func¬ 
tion p,{x,t). In Kress’ early paper Q on the combined field operator, he shows 
fhaf fhe opfimal resulf (for fhe disk) when k is close to zero resulfs from faking 
6 = — I + 0{k‘^) in our nofafion. This is somewhaf af odds wifh our choice to fake 
b real and positive. 

Careful examination shows fhaf Kress’ choice only works if Ai > 0 : In fhe sef 
k 2 < I All fhe real parf of fhe quadratic form in (fTdl) is indefinife (regardless of fhe 
values fhaf a and b fake) and so fo obfain fhe desired resulf we need to employ fhe 
imaginary parf of fhe quadratic form. If 6 = 6 i + ib 2 , then this form would become 



(16) 


In order for the term in the brackets to be definite where A 2 > 0, we need to 
take a > 0. In order for this expression to be definite in both components of the 
set 0 < A2 < I All it is clearly necessary to take 62 = 0 . Of course one could 
use functions of A more complicated than oA + ib as the coefficient of Sk, but 
this would considerably complicate the relationship between the solutions in the 
frequency and the time domains. 

3 Scattering Theory 

Let {a{x),b{x)) be positive functions defined on dD. In fhe recenf paper of Baskin, 
Spence, and Wunsch, see Theorem 1.10 of [T], if is shown fhaf fhere is a positive 
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number /3i so that the operator (A + k‘^)u acting on data in D that satisfies the 
boundary condition {dnU — i{a{x) + ib{x))u) = 0 , is invertible for Im k > -/3i. 
This result does not assume that dD is non-trapping. Using this result, along 
with modern refinements of theorems of Lax and Phillips and Lax, Morawetz and 
Phillips |18]|3, we can prove the following result: 

Theorem 2. For n G N, let D be a non-trapping region in If a and b are 

both positive, then there is a positive number a so that A{a, 6; k) is invertible in 
Irak > —a. 

Proof. The facts we use, in addition to the result of |Tl are 

1. If L) is a non-trapping region then the generator, B, of the compressed wave- 
semigroup Z{t) of Lax and Phillips has its spectrum, ( 7 ( 5 ), in a half plane 
of the form Im k < -/3o < 0 . 

2. A number k belongs to (t{B) if and only if there is an “eventually outgoing” 

solution v+ to (A -I- = 0 with v+ \dD= 0. See Theorem V-4.1 in ||9l. 

By definition (see [911 pp. 126-7) a solution is eventually outgoing (even for Im k < 
0 ) if and only if it can be represented in the form 

v+{x) = / [duyV+{y)gkix - y) - v+{y)d^ygkix - y)]dy, (17) 

dD 

with Uy = —Uy, the outer normal to D'^. See Theorem IV-4.3 in ||9l. Here guix) is 
the “outgoing” fundamental solution for A -|- in 

We argue as before: if for some k, with Im k<0, the operator A{a,b]k) has a 
non-trivial null-space, then so does A'{a, 6 ; k). Let yp 7 ^ 0 be in this latter nullspace 
and set 

tt = y pgkix - y)dSy. (18) 

dD 

Unlike the case where Im A; > 0, we do not know, a priori, that u- 7 ^ 0. Indeed, 
there are just the two cases to consider: u- = 0 or not. If u_ 7 ^ 0, then u- is a 
non-trivial solution to the boundary value problem considered in {T}, and therefore 
Im A: < —/3i. 

The other possibility is that tt_ = 0. This implies that n+ \dD= 0 as well, since 
u is continuous across dD. However need not be zero, since it is exponentially 
growing at infinity. Indeed, if y? 7 ^ 0, and rt_ = 0, then the jump conditions for 
dnU± show that cannot be zero. Because \dD= 0) and dnU+ = —p, 
it follows that n+ satisfies (fTTI) and therefore that n+ is an eventually outgoing 
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solution in the sense of Lax and Phillips. Henee k must belong to the speetrum of 
B, and therefore Im A: < —/3o- This shows that if A{a, b; k) is non-invertible, then 
Im/c < max{—/?!, —/3o} = —a- D 

The eonverse is also true. 

Theorem 3. If k is a scattering pole for the Dirichlet Laplacian, then A{a, 6; k) 
has a non-trivial null-space. 

Proof The theorem of Lax and Phillips states that /c is a seattering pole if and only 
there is an outgoing solution to the BVP 

(A + k^)u = 0 in and u \bD= 0; (19) 

the outgoing eondition in this ease is equivalent to 

u{x) = - J Pkix - y)dnu{y)dSy for x G (20) 

bD 

If we let ip = dnU, then vanishes on bD and 

I - D'kip = p, (21) 

from whieh it is immediate that 

Afa, b] k)p = 0, (22) 

whieh proves the elaim. Thus every seattering resonanee oeeurs among the set 
{k : ker A{a, b; k) 0}, for any ehoiee of positive a and b. □ 

Remark 2. These arguments apply generally to identify the frequeneies for whieh 
the null-spaee of any eombination of single and double layer potentials is non¬ 
trivial. The only eases not explieitly eovered are those of the single and double 
layers alone. One easily establishes that Sk fails to be invertible for in the 
Diriehlet speetrum of D union with the seattering poles of the Diriehlet operator 
in O. The (exterior) double layer Id /2 + Dk fails to be invertible for k'^ in the 
Neumann speetrum of D union with the seattering poles of the Diriehlet operator 
in O. 

Spenee et al. show that that the resolvent kernel for the interior impedanee 
problem satisfies the norm estimate \\Ra^bik‘^)\\ < i+\k\ k 'ma strip around the 
real axis. It is not elear what the analogous bound is for the analytie eontinuation 
of the exterior Diriehlet problem ||ZiDir((i^^)|| into the lower half plane. This oper¬ 
ator maps eompaetly supported data into funetions that have exponential growth at 
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infinity. For later applications we will need to determine bounds on \\A{a,b; k)~^\\ 
for k in the lower half plane. In particular it will be important to establish a bound 
of the form ||^(a, 6; k)~^ || < M(1 + |A:|)™. 

Suppose that g{x,t) is data for the wave equation on the boundary of D with 
support for f G [0, T], and let 

OO 

g{x,k) = j g{x,k)e^^dt. (23) 

0 

The source term on the boundary as a function of t will be given by the inverse 
Fourier transform 

OO 

fi{x,t) = — [ e^^^^^[A{a,b;k)]~^g{x,ik)dk. (24) 

27r J 

— OO 

If this is the case, then the corresponding time-domain representation of a solution 
to the wave-equation is 


u{x, t) = Vtfi + St{adt + b)g, (25) 

with St and Vt the retarded potential single and double layers for the wave equa¬ 
tion. 

Assume that for any e > 0 and Im/c > — (a — e), we have an estimate like 

\\[A{a,b;k)]-^\\<CS + \k\r- 

For g with sufficient smoothness in t we can deform the contour to conclude that, 
for cj < a, we have the representation: 

OO 

g{x,t) = — f e^^^^~^\A{a,b-, k — ia)]~^g{x,ik — a)dk, (26) 


and therefore 

M-,t)\\<Ce-^%\\m,T. (27) 

For conveniently defined norms, this is essentially what we wanted to prove. 

Of course this may not give the “optimal” rate of decay. It should be noted 
that the rate of decay a is directly related to two spectral invariants, one being the 
genuine scattering resonances, and the other the spectrum of the BVP for A on 
LP‘{D) with impedance BC 


dnU — i{ak -b ib)u = 0 on dD. 


(28) 



At least near the negative imaginary axis, we can push this down into the lower half 
plane so that the first eigenvalue encountered is a genuine scattering resonance. 
This follows from a careful examination of the integration by parts argument given 
above. 

Let a and b be positive constants. We assume that (An + k‘^u) = 0, and 

dnU — i{ak + ib)u |"aD= 0. (29) 

If we integrate by parts we see that, with k = ki + ik 2 ■ 


0=1 ((An + k^u)udx 


D 


= — j [\^u\^ + {k2 — ki)\uY]dx — {ak2 + b) J \ufdS+ 

dD 


D 


iki 


j 2 A ;2 / \u\'^dx + a / |np(iS'j . 
\ D dD / 


If ki ^ 0, the we can substitute from the imaginary part the fact that 


/ 

an 


u\^dS = —^ [ InPdx 


D 


into the real part to see that 

0 = / |Vnpdx — 


a{k‘l + k^) + 26^2 


InPdx. 


(30) 


(31) 


(32) 


D D 

Evidently the quantity in the brackets must be positive, which implies that 

kl+(k2 + -) > (- 


(33) 


Thus the impedance boundary value problem has a pole free region lying in the 
disk minus the imaginary axis 

kf + ^2 + —^ < (a) ^ ~ 

Along the imaginary axis its clear that k 2 < — 6/a is needed for (1^ to hold. A 
somewhat better result can be obtained by using the fact that there is an estimate of 
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the form 


/l 

ufdS < Cl 

f [ Vu|^ + Ittpjdx 

J 

dD 


J 

Li) J 

which implies that, if ki 

= 0, then 



(35) 


k 2 < max{r : 1 + (ar + h)Ci < 0 and + Ci{aT + b) < 0}. (36) 

Taking b/a large, we ean arrange to have the pole-free region eneompass as mueh 
of the lower half plane as we want. 

4 The Round Sphere 

In this seetion we eonsider the ease that T is equal to the unit sphere in We 
take a to be a eonstant a and b a eonstant /3. Sinee everything commutes with the 
action of the rotation group, we can analyze this problem one spherical harmonic 
subspace at a time. 

4.1 Integral Equations on Spherical Harmonic Subspaces 

We need to analyze the action of the operator 

, 1 f ny{x-y) f , dty{y,T)\ , 

G - J J dS,+ 

r 

Att J \x-y\ 

r 

where t = t — \x — y\, on data of the form y{y, t) = Y^{y)f{t). Notationally it’s 
easier to study the general class of operators of the form: 

K{Y^f){x,t) = ^ j Yn{y)f{t -\x- y\)k{x ■ y)dSy. (38) 
r 

Up to taking time derivatives of /, all our operators take this form. 

The key observation is that if we let x-*- be a unit vector orthogonal to x, then 
we can define a coordinate system {6, cp) on the unit sphere by setting 

y = cos 9x + sin OR^x^. (39) 
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Here, Rff, is a rotation through angle (p about x. Substituting, we see that 


TT 27r 


KiY:rf)ix,t) = 


Att 


Y^{cos 9x + sin OR^x^ 


0 0 


f{t — a/2(1 — cos9))k{cos 9) sin 9d(j)d9. (40) 


A simple calculation shows that 

2tt 


1 


J Ypp^{cos9x + sm9R,f)X'^)d(j) = Ypp^{x)Pnicos9). (41) 


Here Pniz) is the degree n Legendre polynomials normalized so that P„(±l) = 
(±1)’^. Therefore we have that 


K{YPP'f){x,t) = — / Pn{cos9)f{t— -^2(1 — cos 9))k{cos 9) sm9d9. 


(42) 


Let s = -^2(1 — cos 9), and k{s) = -^2(1 — s)k{s), to obtain: 


i^(yr/)(x,f) = 


y-(x) 


Pn{l - sV2)/(t - s)k{l - s^/2)ds. (43) 


To represent acting on data of this type we need to determine what k is 
for each of the terms in (l37]) . The kernels are 


ny{x- y) ^ x-y-1 - ^ ^ 

[2(l-x.y)]i 2 


ny{x- y) _ x-y-l ^ ^ ^2(1 - s) 


|x-y|2 [2{l-x-y)] 

1 1 


(44) 




- y\ ^2(1 -x-y) 


o- k{s) = 1 
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Acting on Yj^f we get 

2 

/ + 2 y “ 'S^/2){/(i - •s) (-1) + 

0 

dsf(t-s){s)-2adsf(t-s) + 2/3f(t-s)}ds . (45) 
Note that we have replaced dt with —dg. Integrating by parts with respect to s gives 

= iTw 2^/(*) + - 2)- 

2 

i y [(2 - 2/3)Pnil - s^/2) - s{s - 2a)P'(1 - ^2/2)]] f{t - s)ds 
0 

Denote the 1-d operator in the brackets by Gn'^ f- If we let a = /3 = 1, we get a 
very simple integral equation of the second kind to solve for / : 

2 

= fnit) + 11 2)^n(l - sV2))/r(f - s)ds = (47) 

0 

If n = 0, then this is simply /q (f) = 5'o(f). 




4.2 The Roles of a and fd 

The equation 

- 2) = 5(f), (48) 

is quite informative as regards the rate of decay and regularity of the solution. 
Assuming that 0 < a < 1, we simplify this equation to obtain 

1 — O' 2 

fit) + T—fit - 2) = -T—git). (49) 

1 + a 1 + a 

If g{t) = 0 for f < 0, then we can let f{t) = 0, for f < 0 as well. The solution to 
this equation is formally given by the infinite series 


2 1 — rv 

fit) = ^i->^ygit - 2j), with A = 

1 + a 1 + a 

j=0 
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(51) 


For data supported in [0, oo) this sum is finite for any t. In partieular: 

f{t) = for t € (- 00 , 2). 

If g is supported in [0, 2M] then the sum ean be made more explieit: 

2 

fit) = i-^)^9it - 2j). (52) 

i=[f-M]A0 

Eventually / satisfies an esfimafe of fhe form 

\fit)\<CM,x(^\^y Ml^. (53) 

This shows how fhe ehoiee of a affeefs fhe deeay of fhe souree ferm on fhe bound¬ 
ary. Furfhermore, if g is smoofh buf ifs supporf does nof lie in an inferval of fhe 
form [21,2{l + 1)], fhen fhe solution / fypieally has a jump diseonfinuify af every 
posifive even infeger. 

The n = 0, /3 = 0 ease ean be solved by using Laplaee fransform: 

= 7o ('r)^/('r), (54) 

where 

70^ = ^ 1 + (1 - e"^^) . (55) 

Nofe, however, fhaf no matter whaf value a fakes, fhis mulfiplier has a roof af 
s = 0. This explains why we need fo fake /) > 0 in order fo gef exponenfial deeay. 

4.3 Numerical Illustrations 

To provide eonerefe examples of fhe implieafions of our analysis, we have earried 
ouf represenfafive numerieal simulafions in fhe ease n = 0. We eonsider fwo 
ehoiees for fhe Diriehlef dafa: 

g{t) = 8sin(50f) • Oscillatory, 

g{t) = Non — oscillatory. 

We also eonsider three ehoiees for the parameters in the integral equation: 

a = a = 0, b = (3 = 0. 

a = a = 1, b = j3 = 0. 

a = a = 1, b = /3 = -. 
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(Note that the optimal ehoice for the sphere, a = b = 1, leads to a trivial equation 
when n = 0 so we avoid it.) 

Our numerical method is based on the standard Adams predictor-corrector idea. 
Recall that for this special case the equation to be solved is 


/i(f) -f 


1 — a 
1 + a 


/i(f - 2) - 


1-/3 

1 “h Oi 


L 


* 2 
Ks)ds = —— 

t-2 1 + a 


g{t). (56) 


Introducing a time step At we rewrite in correction form: 

1 — a 


lj.{t + At) = /x(t) — 

-f 
-f 


1 “h rr 

1 _^\ rt+^t 

1 “h Oi 

2 


)i: 


(/x(t + At — 2) — /i(t — 2)) 
1-/3^ 


g,{s)ds — 


1 rr 


IJ.{s)ds 


lt-2 


1 cx 


{g{t + At) - g{t)). 


(57) 


In the examples which follow we interpolate past solution data on an interval 
(jAt, (j + l)At) by the Lagrange polynomial of degree 5 using the approximate 
solution at t = kAt, k = j — 3,j + 2. Precisely the interpolant is used to 
calculate p(t — 2), /i(t + At — 2), fj,{s)ds appearing on the right hand 

side of (l57l) . The future integral uses the 6th order Adams-Bashforth-Moulton 
predictor-corrector method. That is, we compute a predicted value + At) 

using the degree 5 Lagrange interpolant of g{t—kAt), A: = 0,..., 5 to approximate 
g{s)ds. We then compute fi{t+At) using the interpolant of (t+At) and 
g{t — kAt), A: = 0,..., 4 inside the integral. We have verified that convergence at 
6th order is generally obtained (in some experiments the observed rate was reduced 
to 5 when a = /3 = 0) and we have also cross-checked the results with those 
computed using 2nd and 4th order solvers constructed the same way. The solutions 
displayed in the graphs were calculated with a time step At = 97/6400 for the 
non-oscillatory data and At = 97/12800 for the oscillatory data. Comparisons 
with a coarsened computation suggest that the solutions are accurate to at least 7 
digits in the non-oscillatory case and 4-6 digits when the data was oscillatory. Here 
the maximum errors in the oscillatory case are two orders of magnitude larger when 
a = /3 = 0; looking at the solutions below we attribute the difference to dispersion 
effects. 

Solutions of the integral equation for the choice a = /3 = 0 are displayed in 
Figure [T] We see that in the case where the data is oscillatory g does not decay but 
oscillates at high frequency. When the data is non-oscillatory we observe linear 
growth. In each case we note that the solution of the wave equation itself decays 
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Figure 1: Density ^{t) at mode 0 with oscillatory data (left) and non-oscillatory 
data (right) for a = /3 = 0. 

exponentially in time, so, as our analysis shows, the long time behavior of the 
density /i(f) is quite different than that of the function it represents. 

Solutions of the integral equation for the choice a = 1, /3 = 0 are displayed in 
Figure |2] We see that in the case where the data is oscillatory ^ apparently decays 
in time. However a closer look indicates that it in fact approaches a very small 
steady state. When the data is non-oscillatory fj, clearly approaches an 0(1) steady 
state, again in contrast with the solution of the wave equation itself. This behavior 
corresponds to the pole at the origin discussed above. 


Solution with osciiiatory data: asl, |ts:0 


■1 





10 "® 10 '’ 10 ° to' 10 ® 


Soiution with nonosciliatory data: Oisl, ^==0 



Figure 2: Density /i(f) at mode 0 with oscillatory data (left) and non-oscillatory 
(right) for a = 1, /? = 0. 

Lastly we consider the solutions which arise when a = 1, /3 = ^; these are 
displayed in Figure [3l Here in each case we have exponential decay of n in time. 
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Solution with nonoscillatory data: a=1,3=1/2 



Figure 3: Density ^{t) at mode 0 with oscillatory data (left) and non-oscillatory 
data (right) for a = 1, /3 = ^. 


5 High Frequency Asymptotics for Convex Bodies 


In the previous section we saw that, at least for round spheres, a principal deter¬ 
minant of the decay rate for the solution of the integral equations on the boundary, 
Gn'^ f = g, is the delay term |[(1 -f a)f{t) + (—1)"'(1 — a)f{t — 2)]. Choosing 
a = 1 removes the sharp contribution from the antipodal point, and gives the opti¬ 
mal rate of decay available for this case. In this section we study the high frequency 
behavior of Dk — i{ka + ib)Sk for a smooth surface F. In general one cannot ex¬ 
pect to have the sort of exact cancellation attained for a round sphere, however we 
demonstrate the rather remarkable fact that for convex bodies, the leading order 
part of the delay term is optimally canceled by taking a = 1. It is less clear how to 
choose 6, but this is related to insuring invertibility of A{a, b] k) at zero frequency. 

We begin by recalling the formula for the relevant kernels: 


Skf{x) 

Dkf{x) 



r 


ik\x-y\ 


^ik\x y\l^(^x — y),ny) 


4:'k\x — yp 


1 

x-y\ 


ik 


f{y)dSy. 


(58) 


As is well known, these are weakly singular integral operators on a surface in 

m3. 

Since we are interested in the high frequency asymptotics, the contribution of 
the diagonal singularity is of less interest to us than the that of the other critical 
points of the phase function 


(px{y) = lx - y\. 


(59) 
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For completeness we state the contribution from the diagonal: 

sty{x) = -M+o{k-^) 

where H{x) = ^(k:i(x) + K2(x)) is the mean curvature of F at x with respect to 
the outer normal vector Uy. Observe that the diagonal asymptotics of the combined 
field operator are given by: 

6; k)fix) = ^f{x) + {2H{x) - + 0{k-^) (61) 

Now we turn to the other critical points of (j)x '■ 


Cx = {y eV ■. Vy(j)x{y) = ±ny} \ {x} 
= {y GT :x = y=f cj)x{y)ny} \ {x}. 


(62) 


Even for a strictly convex body, these critical points can be degenerate. For ex¬ 
ample, if r contains an open subset, U, of a sphere, whose center xq lies on E, 
then U C Cx^ ■ To use the standard techniques of stationary phase we restrict our 
attention to the case of where T is C°° and Cx consists entirely of non-degenerate 
critical points. There are, for example, many boundaries for which Cx consists of 
a global maximum for every x G T. A simple analysis of this property, using the 
characterization of Cx in (l62l) . shows that it is stable under C^ small perturbations. 

If X is a non-degenerate critical point of (j)x, then x lies along the line {x+trix}- 
We choose orthogonal coordinates so that x = 0, nj = (0,0,1), so that x = 
(0,0, id), with d = |x — x|. For a convex body we always have —d. In these 
coordinates we represent the surface F near to 0 as a graph over its tangent plane 
by taking 

y = {z, h{z)) with /i(0) = V/i(0) = 0. (63) 

Fet be a smooth function equal to 1 in a neighborhood of 0, with a single critical 
point of (l>x in its support. 

With X = (0,0, ±d), the asymptotic contribution of x to Skf{x), which we 
denote by S^f{x), is given by 

r „ik{\z\'^^ 2 dh{z)+h{z)^+(f)^ f(z)lb(z)dz 

svix) ~ ^- nzmz)az 

J 47r(|z|2 T 2dh{z) + h{zf + d2)l y^l + |Vd(z)|2 

To apply the stationary phase formula (at least to leading order) we need only 
compute the Hessian of the phase function at 0; it is 


nj :fdhi2i0) \ 

d\Tdh2l{0) lTdh22{0))- 


(65) 


17 








According to stationary phase we see that 


■n-isgn'H g(0) 






2A:Y^det Uxxi^) d 


+ o{k- 


( 66 ) 


Recall that sgn'H3,j(0) is the signature of the Hessian, which is the number of 
positive eigenvalues minus the number of negative eigenvalues. For a maximum 
sgn7f3;j(0) = —2. Even for strictly convex bodies there can be several critical 
points with different signatures. If this happens then there are also points x so that 
Cx contains degenerate critical points. 

The calculation for the double layer is similar: 


mi^) 


^iki\z\^T2dhiz)+h{zr+d^p ^ ^ d)f{z)i^{z) 


dvrd^P =F 2dh{z) + h{zy + 
1 


{\z\^ =F 2dh{z) + h{zY + ^ 2)2 

Once again, applying the stationary phase formula we see that 


I — ik 


dz. (67) 


DUix) 


=Fie*“/(x) 


+ o{k-^). 


( 68 ) 


2Vdet7f,s(0) d 
As noted above, if F is convex, then x = (0,0, —d), and therefore, we see that 


Dlf{x)-ikSlf{x) = 0{k-^), 


(69) 


for all X and non-degenerate critical points x ^ Cx- Unless each Cx consists of a 
single maximum, for every x G F, then, for some points x € T, (j)xi^ sure to have 
degenerate critical points. If these are isolated and of finite degeneracy, then this 
difference is still likely to be o(l) as /c —> oo. Of course if F is not convex, then it 
is perfectly possible that we may sometimes get the — sign in the asymptotics for 
Df.f{x), equation (IMI) . so that we obtain contributions of the form: 

Dlf{x)-iakSlf{x) =~ -i(l + a)^====— -^ + 0(fc-i). (70) 

Vdet7fa;j(0) |x-x| 

Since a > 0 is needed to keep poles out of the upper half plane, this may prevent 
any source for the wave equation, defined by solving dUl wifh an allowable choice 
of (a, 6), from decaying exponentially. 
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